Resonance fluorescence in ultrafast and intense x-ray free-electron-laser pulses 
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The spectrum of resonance fluorescence is calculated for a two-level system excited by an intense, 
ultrashort x-ray pulse made available for instance by free-electron lasers such as the Linac Coherent 
Light Source. We allow for inner-shell hole decay widths and destruction of the system by further 
photoionization. This two-level description is employed to model neon cations strongly driven by 
x rays tuned to the ls2p _1 — ► Is 2p transition at 848 eV; the x rays induce Rabi oscillations 
which are so fast that they compete with Ne ls-hole decay. We predict resonance fluorescence 
spectra for two different scenarios: first, chaotic pulses based on the self- amplified spontaneous 
emission principle, like those presently generated at x-ray free-electron-laser facilities and, second, 
Gaussian pulses which will become available in the foreseeable future with self-seeding techniques. 
As an example of the exciting opportunities derived from the use of seeding methods, we predict, in 
spite of above obstacles, the possibility to distinguish at x-ray frequencies a clear signature of Rabi 
flopping in the spectrum of resonance fluorescence. 

PACS numbers: 32.50.+d, 32.30.Rj, 42.50.Ct, 32.70.Jz 
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I. INTRODUCTION 

The spectrum of resonance fluorescence which is emit- 
ted by an ensemble of atoms and ions driven by an intense 
near- resonant electric field [1, 2] is one of the cornerstones 
of quantum optics [3, 4]. The spectrum is measured ex- 
perimentally by exposing an atomic ensemble to intense 
light and detecting the scattered photons as shown in 
Fig. 1. During the last few decades such studies have 
received wide attention and have stimulated the develop- 
ment of non-perturbative methods in quantum electrody- 
namics for the study of the coherent interaction between 
light and matter [5-10]. 

The resonance fluorescence spectrum of a two-level 
system driven by a monochromatic electric field is the 
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Figure f. (Color online) An atom ensemble (red) is driven 
by x rays (blue) tuned to a resonance. The emitted photons 
(green) are measured perpendicularly to the propagation di- 
rection of the x rays. 
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simplest case and has been studied extensively at opti- 
cal frequencies [5, 6, 11-13]. For a sufficiently strong 
continuous- wave (cw) driving field a nonlinear three-peak 
structure appears in the spectrum [14-16] which is ex- 
plained theoretically by the nonperturbative approach of 
Mollow [5, 6]. The presence of three peaks, frequently 
called dynamic Stark splitting, is explained as the result 
of the dressing of bare levels by the external field [17]. 

A cw field is one of the few cases for which an exact an- 
alytical solution of the equations of motion (EOMs) of the 
two-level system exists. When the system interacts with 
a short pulse, a special class of time-dependent functions, 
including the case of a hyperbolic secant pulse, were an- 
alytically explored for particular values of the physical 
parameters [18-20] and a rich multipeak structure in the 
spectrum of resonance fluorescence was predicted [21-27]. 
This property, which still represents a signature of Rabi 
oscillations induced by the intense driving field, is also 
predicted to depend upon the pulse area, but cannot be 
intuitively explained by means of dressed states [28] . 

In this paper we investigate, in terms of a two-level 
model, the coherent interaction of x rays with core elec- 
trons by exciting AT-shell transitions. Previous studies 
of strong-field resonance fluorescence have been relevant 
only at optical frequencies, for which a wide range of 
models and schemes have been investigated [29-35], be- 
cause of the lack of coherent and sufficiently intense light 
sources at short wavelengths. The recent construction of 
x-ray free-electron lasers (XFELs) [36-40] provides one 
with tunable x-ray pulses of unprecedented brilliance, up 
to one billion times higher than the intensity available at 
third-generation synchrotron facilities. The intense and 
ultrafast pulses now available at XFELs offer the oppor- 
tunity to study nonlinear physics at short wavelengths 
[41-50]. In the particular case that we are going to in- 
vestigate here, x rays are able to induce stimulated emis- 
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sion and absorption of photons (Rabi flopping) at a time 
scale that can be compared to and, therefore, compete 
with the ultrafast inner-shell Auger decay [51, 52]. 

Existing facilities such as the Linac Coherent Light 
Source (LCLS) are based on the principle of self-amplified 
spontaneous emission (SASE) [53, 54], i.e., the beam shot 
noise gives rise to the emitted radiation which, as a result, 
possesses only partial temporal coherence and a spiky 
temporal profile. An analogous situation occurred at the 
beginning of optical laser science, when it was timely to 
study the interaction with the chaotic pulses available 
at that time [55, 56]. Self-seeding or optical laser-seeding 
methods are being developed, for which the emitted light 
is produced by the amplification of a regular (Gaussian) 
seeding pulse, which exhibits high temporal coherence 
[57-59]. The rapid development of XFEL sources makes, 
therefore, further theoretical work timely [60]. 

In a recent experiment, intense and ultrashort x-ray 
pulses from the LCLS have been used to excite the 
Is 2p- x -> Is -1 2p transition at 848 eV in Ne + [52]. The 
electron spectrum of resonant Auger decay was measured 
to investigate Rabi flopping. With only partial coher- 
ence of the SASE pulses presently available at LCLS and 
the lack of means for single-shot diagnostics, though, the 
clear observation of Rabi oscillations and its distinction 
from noise effects is challenging [52]. 

Auger decay is the predominant mechanism of inner- 
shell hole decay of light elements such as neon. Because 
of the large Auger yield, high-resolution electron spec- 
troscopy is well suited for soft x rays. Besides Auger 
decay, K holes also decay by x-ray fluorescence, i.e., by 
spontaneous emission of photons, while the system is 
driven by an external field. The spectrum of resonance 
fluorescence represents an alternative way to study the 
coherent and nonlinear interaction between x rays and 
atoms and ions. It complements the results coming from 
the detection of electron spectra of resonant Auger decay. 
High gas densities can generally be used for x-ray emis- 
sion spectroscopy — orders of magnitude larger than with 
electron spectroscopy — which compensates for the small 
fluorescence yield and enables high-resolution measure- 
ments with gratings or crystal spectrometers. In the case 
of resonance fluorescence, however, self-absorption can 
produce line broadening [61], so the gas density and path 
length need to be adjusted to minimize self-absorption 
effects and make them negligible. Furthermore, photons 
are scattered much less off electrons or ions in the in- 
teraction volume than electrons, i.e., space-charge effects 
are of little concern [62] . 

For x-ray energies, present instruments are expected to 
detect the fluorescence spectrum with high frequency res- 
olution. For the purposes of this paper, there are at least 
three choices of instruments: a cryogenic spectrometer 
[63] , a diffraction grating [64] and a crystal spectrometer 
[65]. Cryogenic spectrometers such as microcalorimeters 
are mostly used in the astrophysics community for de- 
tecting x rays from atoms with high atomic numbers and 
high fluorescence yield [66]; at 1 keV, Ref. [67] suggests 



a frequency resolution lower than 0.8 eV. For modern 
grating instruments based on the design described in [68] 
and [69] a resolution of 0.4 eV at 848 eV is expected. It 
might also be possible to achieve higher frequency resolu- 
tion by using higher-order reflections from gratings and 
crystals — though with loss of detection efficiency. With 
the use of wavelength-dispersive spectrometers, such as 
diffraction gratings and crystal spectrometers, one might 
take advantage of their polarization sensitivity in a par- 
allel and perpendicular setup, e.g., for background reduc- 
tion [70]. 

In this paper, we develop a time-dependent theory of 
resonance fluorescence to study the interaction of a two- 
level model with x-ray pulses. In Sec. II we describe our 
theoretical approach, by defining the spectrum of reso- 
nance fluorescence and its main properties and by intro- 
ducing the two-level model that is used throughout the 
paper. Results are discussed in Sec. Ill, where the spec- 
trum of resonance fluorescence is examined for different 
driving pulses. In particular, we compare different spec- 
tra for the presently available chaotic pulses produced via 
the SASE principle and for pulses with a Gaussian tem- 
poral profile that seeding techniques are making avail- 
able. Section IV concludes the paper. Atomic units are 
used throughout unless otherwise stated. 

II. THEORETICAL MODEL 
A. Two-level model 

The coherent interaction between atoms and ions and 
x rays tuned to a particular atomic resonance can be de- 
scribed in terms of a two-level model when the transition 
is isolated from other levels. In our case, we use such a 
model to study the Is 2p~ x — > Is" 1 2p z transition in Ne + 
at an energy of 848 eV [52], driven by a near-resonant 
electric field linearly polarized along the z direction. The 
two-level model, which is depicted in Fig. 2, is justified 
by the fact that the transition is very well isolated, by 
more than 70 natural linewidths separated from the next 
Rydberg excitation Is — > 3p of neutral Ne at 867 eV [52]. 
For neon, relativistic effects and fine-structure splitting 
do not play an important role and, therefore, spin-orbit 
splitting can be neglected. 

We describe the emitted fluorescent light field by a 
quantum operator E(r, t) = E + (r, t) + E~(r, t), where 
E + (r, t) and E~(r, t) are respectively the positive- 
frequency and negative-frequency parts of the operator 
[71]. However, it is sufficient to describe the relatively 
strong driving field classically [12], via 

£{t) = £ (t) cosKt + tfix(t) + ¥>x,o], (1) 

where £o (t) is the time-dependent electric- field envelope, 
ljx is the x-ray central frequency, fx(t) is the time- 
dependent phase of the field, and (px,o 1S the carrier- 
envelope phase (CEP). We assume throughout an elec- 
tric field linearly polarized along the z direction, £o(t) = 
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Figure 2. (Color online) Two-level model used to describe the 
coherent interaction between Ne + and the external driving 
field tuned to the Is 2p _1 — > Is -1 2p transition at 848 eV [52]. 
The ground state Is 2p~ is given by |1) = \L = 1, Ml = 0} 
and the core-excited state Is" 1 2p is written as \2) = \L = 
0, Ml = 0). The external field is linearly polarized along the 
z direction and induces Rabi flopping between these states. 
Spontaneous decay, however, also allows the core-excited state 
to decay to valence-ionized states with Ml — ±1. 



£q CO &z ■, with the unit vector in z direction e z . The 
use of planar undulators at XFEL facilities, in fact, pro- 
duces linearly polarized x-ray pulses [40]; experimental 
evidence for a very high degree of linear polarization of 
LCLS x rays has been given in Refs. [42, 52]. We fur- 
ther assume a pulse with uniform intensity distribution 
profile; spatial averaging is therefore not performed. 

In order to properly model the atomic transition, we 
see in Fig. 2 that the Is 2p _1 configuration is a spin dou- 
blet state with a total orbital angular momentum of L = 
1; consequently, it is triply degenerate in energy. The 
three eigenstates of the unperturbed atomic Hamiltonian 
Hq with energy u>\ , which diagonalize the z component of 
the total orbital angular momentum operator, are |1+), 
| 1q) and |1-), respectively with M h = +1, 0, — 1. Con- 
versely, the Is -1 2p configuration corresponds to the sin- 
gle eigenstate |2) of the field- free atomic Hamiltonian, 
with L = 0, Ml = 0, and energy ui2- The energy of the 
atomic transition is u>2i — U2 — uji- The relevant raising 
and lowering atomic operators are 



interaction of the ion with the classical, linearly polar- 
ized, near-resonant field (1) [3], 

flint = -P ■ £ (t) cos[u x t + ipx(t) + <£ X ,o]- (4) 

The operator P in (4) is the total atomic polarization 
operator 



P = P+ +P 
with P = [P+]t and [72] 



(5) 



1+|P|2) a 1+2 + (1 |P|2) <7i o2 + (1-|P|2) r>i_2 



p(e a + (Ti +2 



(6) 



where e x , e y , and e z are unit vectors in the x, y, and 



directions and e„± 



2 are circular po- 



larization vectors, positive (negative) for polarizations 
A = ±1, with e a ± = —e*^. Due to spherical symme- 
try, the dipole matrix element is real, p = p* , and is 
the same for all transitions; as the atomic states have 
a definite parity, (l±,o|-P|l±,o) = (2|P|2) = 0, dipole 
transitions only couple states with different total angu- 
lar momentum quantum numbers, AL = 1 [72]. 

Since the external electric field is assumed to be lin- 
early polarized [42, 52], the dipole interaction only cou- 
ples the states |2) and |1 ) satisfying the condition 
AMl = [72]: in Eq. (4), within the rotating- wave ap- 
proximation [3] and by using Eq. (6), Hi nt reduces to 



Hi 



int. 



(a lo2 e i ^^( ( ']+a 21l e -iNtt+*xMl) ; 

(7) 

where we have set the CEP tfx,o to and where the 
instantaneous Rabi frequency 



n R {t) = fp£ (t) 



(8) 



has been introduced. 

In our model the dynamics of the two x-ray coupled 
states 1 2) and | lo) is independent from the other two 
states |1±) and one can develop an actual two-level de- 
scription of the system in which the EOMs exclusively 
contain the two aforementioned states | lo) = |1) an d |2) 
and neglect the other two states entirely. 



% = H)(il, *, 3 e {1+, lo, 1-, 2}. (2) 

The interaction between the ions and the electric field 
is described in the dipole approximation because the Is 
orbital of neon is very compact, involving dimensions 
much smaller than the wavelength associated with the 
transition ls2p _1 — > Is" 1 2p, such that nondipole terms 
are small [51]. The Hamiltonian of the system is 

fl = An + flint , (3) 

where flo = w * ^« ^ s the unperturbed atomic Hamil- 
tonian with eigenvalues u)i, whereas fli nt represents the 



B. Density matrix formulation and equations of 
motion 

We investigate in the following the two-level system 
formed by the states |1) = | lo) an d |2). 
We introduce the density matrix 

Pij(t) = = (%( f )> (9) 

(i, j G {1, 2}), whose evolution is described by the mas- 
ter equation [73] 

% = -\[H,p{t)}+cp(t)+mt)- (io) 
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The first term — i [H, p(t)] describes the coherent dy- 
namics of the two-level system. In the total Hamilto- 
nian H [Eq. (3)] the only relevant terms of the unper- 
turbed atomic Hamiltonian H Q are u>i<7n and uj 2 ^22 ■ The 
Lindblad operator Cp(t) represents the norm-conserving 
spontaneous decay of the population from the excited 
state 1 2) to the ground state The rate at which this 
process occurs is given by [3] 



4u^ 



3c 3 



(11) 



Atoms and ions with high atomic numbers are usually 
characterized by a high fluorescence yield, i.e., the im- 
portance of spontaneous decay increases with the atomic 
number of the ion of interest. The last term T>p(t) de- 
notes the norm-nonconserving term not present in the 
Lindblad form of the master equation [74]. We intro- 
duce this term to describe the decrease of the population 
of both the upper and lower states [4]. These norm- 
nonconserving processes include Auger decay, photoion- 
ization of the system due to the intense external field, and 
spontaneous decay from the excited level |2) to the two 
levels which are not coupled by dipole interaction. 
We do not include Doppler broadening [61] and collision 
effects [4] in our model, because they involve time scales 
much longer than the decay time of the system and, at 
room temperature and for a pressure of 1 atm, they can 
be neglected [75]. 

Auger decay and photoionization destroy the two-level 
system by further ionization of Ne + to levels which need 
not be taken into account explicitly. The Auger decay 
width is 2~a, whereas the rate of photoionization Jp i(t), 
1 e {1, 2}, is [76] 



(12) 



with the photoionization cross section for the level i 
0Xj = (Ti{cJx), the x-ray flux 



Jx{t) = J(t)/wx 
and the x-ray intensity 



i(t) 



8na 



(13) 



(14) 



Notice that we evaluate the photoionization cross section 
and the flux at ujx since the cross sections do not vary 
much within the bandwidth of the field. 



The spontaneous decay of the excited level |2) to the 
states \1±) also represents a process which does not con- 
serve the norm of our two-level system. The total radia- 
tive decay rate is given by 



r R = r, 



r, 



R.cr- 



(15) 



where Ir z is the spontaneous decay width to the state 
|1) given in (11) and /r jCT ± are defined analogously for 
the other two decay channels; the second equality exploits 
Eq. (6) . Since the spontaneous decay of the excited level 
1 2) only depends on the population of the state itself, as 
we are going to show in the following EOMs, the actual 
dynamics of states |1±) can be indeed neglected for our 
purposes. 

The total decay processes are included in Eq. (10). 
In order to derive the EOMs for the four relevant com- 
ponents of the density matrix, we move to the rotating 
frame [77], by introducing the operators 



(16) 



(17) 

Pii(t), 



T ii, ?12 — e <Ti2, C21 — e CT21, 



whose expectation values are denoted by 

Rij(t) = (&(i)>, 

which, from (9) and (16), implies that Ru(t) 
R12 = P12 c~ iux * and R 21 = P21 e iux *. 
We introduce the vector 

R(t) = (R n (t), R 12 {t), R 21 (t), R 2 2{t)f 



whose components are given by the elements of the den- 
sity matrix Rij(t) in the rotating frame. Before the ar- 
rival of the light pulse the two-level system is assumed 
to be in the ground state, i.e., Rq = (1, 0, 0, 0) T . This 
assumption is supported by experimental observations of 
orbital alignment in ions produced by strong-field ioniza- 
tion [78]. If the fraction of ions in the Ml = ground 
state is lower than 1, the resonance fluorescence spectrum 
must be multiplied by this factor. 

The master equation (10) can be rewritten in matrix 
form 



dR{t) 
dt 



= M(t)R(t), R(0) = Ro, (18) 



where M(i) is the following time-dependent matrix 



M(t) 



-7i (*) 



V 



•ivx(t) 



1 2 e 



2(7lW 



e -i<px(t) 

72 W) 











-i^-|(7i(*) 



72W) 



j£M£) e i<px(t) 
[ f2 R(*) e -i(Px(t) 

-72 (*) / 



(19) 



with 



7i(*) =vx,iJx(t), 
72(t)=a x , 2 Jx(t)+r A + r R 



(20a) 
(20b) 



where we have defined the detuning A — lj 2 i — ux. 
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The knowledge of the time evolution of the atomic one- 
time expectation values is used to derive the two-time 
expectation values necessary for the computation of the 
spectrum of resonance fluorescence. For this purpose, we 
introduce the two-time vector 

Y(t u t 2 ) 

= (Y n (ti,t 2 ), Y 12 (ti,t 2 ), Y 21 (t 1 ,t 2 ), Y 2 2(ti,t 2 )) T , 

(21) 

whose elements are defined as 

Yijituh) = (^(fiKiate)). (22) 
Applying the quantum regression theorem [77, 79] yields 

dY{t u t 2 ) 



M(ti)r(*i,* 2 ), h>t 2 , (23) 



with the initial conditions given by Yij(t 2 ,t 2 ) = 
SnR2j(t2)- The solution of the first set of differential 
equations (18) provides one with the initial conditions 
for the second set of differential equations (23), whose 
solution gives 



Y 12 (h, t 2 ) = (Ql(*l)«fl2(*2)) 



(«T2l(il)<7l2(i2))e- iWx ( tl - t2 ). 



(24) 



Here, S(ia, r)dtodA represents the energy detected on 
average in the differential energy interval [u>, co + dco] and 
in a surface element dA = r 2 dfi e r centered on r = 
re r . Further, a is the fine-structure constant, dfi is the 
differential solid angle, and e r = r/\r\ is the unit vector 
in the direction of observation (0 is the position of the 
atom). 

We assume that the driving field propagates along 
the y axis. In the far-field limit and in the electric- 
dipole approximation — away from the y propagation axis 
in which also the driving field would be present — the 
electric-field operator associated with the fluorescent 
light can be related to the atomic polarization operator 
P+(t) [Eq. (6)] via the relation [3, 7, 12] 



2 

E+(r, t) = ^f- {P+{t - r/c) - [P+(t - r/c) • e r ] e r } . 

(27) 

If the detector is placed along the x axis, as shown in 
Fig. 1, then e r = e x and one obtains from (6) 

E+ (r e x ,t) = E+ (r e x , t) e z + E+ (r e x ,t)e y , (28) 
with 

E+(re x ,t) = ^a lo2 (t-r/c) (29) 
c r 

and 



C. Spectrum of resonance fluorescence 

The study and computation of the spectral properties 
of the fluorescent light requires the knowledge of the first- 
order autocorrelation function of the electric-field opera- 
tor [80, 81] 

G (1) (ti,i 2 ,r) = (E-{r, ti) ■ E+{r, t 2 )>. (25) 

In the case of cw light, when the first-order autocorre- 
lation function depends explicitly on the time difference 
t = t x - t 2 , i.e., G^(h,t 2 ,r) = G^{r,r), the Wiener- 
Khintchine theorem [3] states that the power spectrum 
of resonance fluorescence associated with the rate of pho- 
tons emitted at a given frequency is well defined and 
given by the Fourier transform of (jW(t, r) [82]. How- 
ever, for ultrashort light pulses, (t±, t 2 , r) explicitly 
depends upon the two distinct instants t\ and t 2 and 
the Wiener-Khintchine theorem cannot be analogously 
used to define a power spectrum. Instead, one needs to 
study the energy spectrum of resonance fluorescence, de- 
fined as a quantity proportional to the probability that 
an ideal photon detector — modeled itself as a two-level 
system with tunable transition energy uj — is excited by 
the fluorescent light. In first order of perturbation the- 
ory and in the electric-dipole approximation, the energy 
spectrum is defined as [80] 



E+(re x , t) = ^= [*i+a(* ~ r / c ) + *i-a(* - r/c)] , 

(30) 

whereas, because of the placement of the detector, the x 
component of the electric-field operator E£ (r e x , t) van- 
ishes. Analogously, the autocorrelation function is split 
into two parts: 

G^(h,t 2 ,re x )=(E-{re x , h)E+(re x , t 2 )) 



+ 00 Z' + OO 



oo J — oc 



G {1) (t 1 ,t 2 ,r)e' iuj( - t ^ t2 '>dt 1 dt2. 

(26) 



+ (E-(re x , t x )E+(re x , t 2 )). 



(31) 



The first term in (31) is the autocorrelation function of 
the fluorescence photons which are polarized along the z 
direction; the transition with which they are associated 
(|1) — » 1 2)) is driven by the external field, which modu- 
lates the polarization operator along the z direction and 
induces Rabi flopping. The general case of a detector 
placed in the x — z plane, forming an arbitrary angle 9 
with respect to the x axis, is discussed in Appendix A. 

With (26) and (31), the resonance fluorescence en- 
ergy spectrum is also split into two terms S{uj, r e x ) = 
S z {uj, r e x ) + S y {uj 1 re x ). The calculation of S y (uj, r e x ) 
goes beyond the two-level approximation we adopt in this 
paper and requires a complete four-level description of 
the system. In this paper we describe the appearance 
of Rabi flopping in the resonance fluorescence spectrum 
for those photons which are emitted in the transition to 
the ground state |1). For e r = e x this represents the only 
contribution in S z {uj, r e x ) and its observation can be ex- 
perimentally realized using a polarization-dependent de- 
tection to selectively detect the radiation which is linearly 
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polarized in the z direction, in order to select those flu- 
orescence photons associated with the transition to the 
state with Ml = 0. 

Polarization-dependent measurements can be very in- 
formative, e.g., they have played an important role for 
molecules where the valence orbitals can be resolved 
[83, 84]. Reflections from mirrors, gratings, or crystals 
at angles that achieve high polarization selectivity at the 
frequency of the atomic transition involved allow one to 
measure the polarization of the radiation. The use of 
wavelength-dispersive spectrometers, which involve re- 
flecting x rays from a grating or crystal, can provide 
one with polarization selectivity. Energy-dispersive spec- 
trometers, such as a cryogenic spectrometer, can be po- 
larization sensitive if they are pixelated and the x rays are 
hard enough to Compton scatter in the absorber. In ad- 
dition, polarization-dependent detection in a parallel and 
perpendicular setup facilitates background reduction. 

By exploiting the polarization properties of the emitted 
light, we can focus on the first component of the first- 
order autocorrelation function (31), which is expressed 
as 



G« (t u t 2 ,re x )= I(r) (a 21 (ij - r/c) a 12 (f a - r/c)) , 



where 



*M " (^) 



(32) 



(33) 



is a factor dependent on the position of observation at 
which the detector is placed and having the dimension of 
an intensity [12]. 

By introducing the time delay r = t\ — t 2 and notic- 
ing that (o- 2 i (£1)0-12(^2)) = (o- 2 i(i2)o-i2(fi))*, we con- 
clude that knowledge of {a 2 \{t\)cr V2 {t 2 )) in the region 
t\ > t 2 (and hence r > 0) is sufficient for the calcu- 
lation of the energy spectrum of resonance fluorescence 
[77]. We rewrite (26) in compact form as 



Sz(u, re x ) 



8?rr 2 j_ OQ jq 
x <7i 2 (i2))]dr dt 2 , 



/+00 /-t-00 
/ Re[e- iuT {a 21 (t 2 +T) 
-00 Jo 

(34) 



where we use Eqs. (11), 
can use Yi 2 (ii, £2) from 
the energy spectrum of r 

In the following, 
S z (u, Q) = r 2 S z (ui,r 
x axis. S z {u>, J?) df2 dco 
and do>; in atomic units 
1/sr. Finally, we notice 
emitted into dfl is 



(32) and (33). As a result, one 
the solution of (23) to compute 
esonance fluorescence, 
we are going to compute 
e x ) for a detector along the 
is the energy emitted into dfi 
S z (uj, f2) has the dimension of 
that the total detected energy 



« = l +0 ° S z (lu, Q) du = 3/ \ z ^ 21 [ +0 ° R 22 (t) dt, 



exploiting the relation 

/+00 
e -iu>(t!-t 2 ) 
-OO 



III. RESULTS AND DISCUSSION 

Here we apply our two-level model to study neon 
cations on the Ls2p~ 1 — > ls~ 1 2p transition at oj 2 i = 
848 cV [52], i.e., the detuning is A = lu 2 \ — ujx = 0. 
Scattered x rays could be observed if the XFEL beam 
energy is detuned from resonance. As demonstrated in 
Ref. [85], for example, Compton scattering, resonant Ra- 
man scattering, and Rayleigh scattering can be observed 
as the resonance is approached from below. At 848 eV, 
however, resonance fluorescence will dominate the mea- 
sured spectrum. 

The destruction rate of our effective two-level system is 
dominated by the Auger decay width of Ne Is" 1 which is 
T A = 0.27 cV [86]. The dipole moment p = 0.0524 a is 
computed with the Hartree-Fock-Slater mean-field model 
[87-89], whereas the photoionization cross sections are 
computed using the Los Alamos Atomic Physics Codes 
[90, 91]. From Eq. (11) the radiative decay width follows, 
where = 0.0012 eV and the total decay rate is Jr = 
0.0039 eV [92], in good agreement with Eq. (15). 

The spectrum S z (lo, ft) that we will compute repre- 
sents the emitted photons linearly polarized along the 
z direction from our two-level model. Off-resonant 
Rayleigh scattering from 2s and 2p electrons in Ne is, 
however, not taken into account. This elastic scattering 
is predicted to be anisotropically distributed for a lin- 
early polarized electric field: in our case, S(t) = £(t) e z , 
the intensity of the elastic scattering would be affected 
by the source-dependent polarization factor sin 2 ip [76] , 
where ip is the angle between the z axis and the direc- 
tion of detection e r at which the detector is placed. This 
additional contribution is not included in the two-level 
approximation that we implement in this paper. Its only 
effect is an enlargement of the central peak of the spec- 
trum. 



A. Gaussian x-ray pulses 

Self-seeding techniques at LCLS are providing one with 
pulses with an approximately Gaussian temporal profile 
[58, 59]; it is interesting therefore to predict the evolution 
of the atomic properties in time and the spectrum of 
resonance fluorescence for Ne + cations for this case. We 
write the Gaussian pulse as 

,-[i 2 /(2T 2 )] 



£o,g(0 — £max e 



= 0, (36) 



where T = r G /(2v1n2) and r G is the FWHM of ££ G (t). 
The FWHM of |£ ,gM| 2 is Aw G = 41n2/r G , where 



(35) 
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is the Fourier transform of £o,g(^) [71]. The peak inten- 
sity [Eq. (14)] is Iq = E^x/ (8ira), yielding a maximum 
Rabi frequency [Eq. (8)] J?R G , m ax = p£ ma x = p\/8iTaI G . 
Further, we introduce the pulse area 

/+oo 
Q R (t) dt, (37) 
-oo 

which was shown to play an important role in the descrip- 
tion of the dynamics of a two-level system in interaction 
with a regular pulse [ipx(t) = 0] and in the properties of 
the corresponding resonance fluorescence spectrum [21- 
24]. Let us assume for now that level decay and pho- 
toionization are both negligible. Then, for A = 0, if 
n is a natural number and Q = 27rn, the final popula- 
tion after the interaction with the pulse is in the ground 
state, whereas for Q = 2-k(ti + 1/2) a complete inversion 
happens and the total final population occupies the ex- 
cited state. For a Gaussian regular pulse the area (37) is 

Qq = ^RG.max Tq y/ 7T / (2 In 2) . 

We begin by studying the interaction of Ne + cations 
driven by a Gaussian x-ray pulse with peak intensity 
Iq — 2.6 x 10 17 W/cm 2 and tq = 5 fs: such x-ray pulses 
will be available in the foreseeable future from seeding 
techniques implemented at LCLS. In Fig. 3 we show the 
time evolution of the two-level system when Auger decay 
is included and when it is not included (Pa = 0) in the 
EOMs (18). The time evolution of the total population 
of the system reveals that Auger decay is the major de- 
population mechanism. Photoionization makes, however, 
also a noticeable contribution at the chosen x-ray inten- 
sity; the maximum rates of photoionization are [Eq. (12)] 
Ppi.max = 0.03 eV and Pp2, ma x = 0.04 eV, which is small 
compared with the Auger decay width Pa = 0.27 eV. In 
Ref. [51] this channel was, therefore, neglected entirely. 
The decay time associated with Auger decay is approxi- 
mately given by At = 1/Pa = 2.4 fs. As one notices in 
Fig. 3, the total population of the system almost com- 
pletely vanishes after the pulse of 5 fs. Whether Auger 
decay is included or not does not interfere with Rabi 
oscillations which are clearly discernible; the pulse area 
Qq = 7 x 2ir results in seven oscillations. 

The corresponding energy spectra of resonance flu- 
orescence are shown in Fig. 4. The Rabi oscillations 
induced by the intense external x-ray field appear in 
both cases with and without Auger decay, with nonva- 
nishing contributions in the region approximately given 
by [-/?RG,max, ^RG.max], with the maximum Rabi fre- 
quency -f^RCmax = 3.9 eV. First, when only spontaneous 
decay and photoionization are taken into account, a mul- 
tipeak structure is predicted, in analogy to what was 
computed in the absence of any decay process [21]. The 
presence of many peaks is nontrivially related to the pulse 
shape of the electric field, i.e., to its finite duration and 
width. The seven peaks in the energy spectrum — six lat- 
eral peaks and the seventh central one — are related, as 
was shown in [21], to the pulse area Qq — 7x2w. Second, 
when Auger decay is taken into account, the multipeak 
structure of the spectrum becomes smoother because of 
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Figure 3. (Color online) Time evolution of the population of a 
two-level system driven by a Gaussian x-ray pulse (36) of peak 
intensity I G = 2.6 x 10 17 W/cm 2 and FWHM duration r G = 
5 fs, corresponding to a pulse area of Qq = 147r. The dashed 
lines show the total population of the two-level system pn(t) + 
p22(t) [Eq. (9)] in the absence (black line) and presence (red 
line) of Auger decay. The solid lines show the corresponding 
occupation of the excited state p22(t). 
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Figure 4. (Color online) Energy spectrum of resonance fluo- 
rescence for the Gaussian pulse used in Fig. 3 in the absence 
(black, dashed line) and presence (red, solid line) of Auger de- 
cay. As indicated by the arrows, the scale on the left refers to 
the black, dashed curve, whereas the scale on the right refers 
to the red, solid curve. 



the increase in the decay rate. Furthermore, the inten- 
sity of the radiation emitted by the two-level system de- 
creases, since Auger decay destroys it and, consequently, 
reduces the fraction of atoms which can Rabi flop. The 
resulting maximum Rabi frequency i^RG^max = 3.9 eV is 
however much higher than the bandwidth of the pulse, 
Auiq = 0.36 eV, the Auger decay width, P A = 0.27 eV, 
and the frequency resolution of present spectrometers, 
Aw rcs = 0.4 eV [68]. Hence the signature of Rabi flop- 
ping, clearly visible in Fig. 4, will be detectable. 

In Fig. 5 and 6 we consider different pulses with 
tq = 2 fs and Qq = 27r(?i — 1/2) [panels (a)] or Qq = 2irn 
[panels (b)], for n 6 {1, 2, 3}. One can clearly see a de- 
pendence of the population of the two-level system upon 
the area Qq. When this area is an odd multiple of tt 
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Figure 5. (Color online) Time evolution of the population 
of the excited state pi2(f) for a two- level system driven by 
Gaussian x-ray pulses [Eq. (36)] of different peak intensities 
(shown in the legend) and a FWHM duration tq = 2 fs. In 
panel (a) pulse areas Qg = 27r(n— 1/2), for n = 1 (red, dotted 
line), n = 2 (black, dashed line) and n = 3 (green, solid line) 
are used. In panel (6) pulse areas Qg = 2nn, for n = 1 (red, 
dotted line), n = 2 (black, dashed line) and n = 3 (green, 
solid line) are used. 
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Figure 6. (Color online) Spectrum of resonance fluorescence of 
a two-level system driven by Gaussian x-ray pulses of different 
peak intensities (shown in the legend) and a FWHM duration 
tg = 2 fs. Line styles of panels (a) and (b) as in Fig. 5. 



[Fig. 5a], a major part of the population at the end of 
the pulse occupies the excited state: one can discern the 
n — 1/2 oscillations due to the interaction with the pulse 
and the following Auger decay of the system when the 
pulse is over. As shown in Fig. 6a, the long Auger decay 
which follows the interaction with the pulse results in a 
high Lorentzian peak in the spectrum of resonance fluo- 
rescence at oj = with a width that can be related to 
the major decay process, i.e., J\. This peak results from 
the fact that the excited system decays freely, radiatively 
and electronically, without Rabi flopping. In contrast, a 
considerably different situation appears when the area of 
the pulse is an even multiple of it [Fig. 56]: after n com- 
plete oscillations, the population of the excited state is 
almost at the end of the pulse. Consequently, the pre- 
viously present post-x-ray-exposure decay does not take 
place. The total emitted energy is therefore lower be- 
cause the central peak at W21 is reduced by almost one 
order of magnitude, as one can clearly see by looking at 
Fig. 66. In the case of a longer pulse, so that the two- level 
system has completely Auger decayed before its conclu- 
sion, the difference between pulses whose area is an odd 
or even multiple of it becomes less important. 

The dependence of the resonance fluorescence spec- 
trum upon the duration of the pulse is an additional 
point that needs to be investigated. In order to observe 
this dependence, in Fig. 7 we study the main features of 
the spectrum as functions of the normalized pulse area 
Qg/(2tt) and of the pulse FWHM duration tq. We re- 
call that for fixed tq the area of the Gaussian pulse is 
directly propor tional to the square root of the intensity, 
Qg = ZftpTG \J {ot/ In 2) Iq. In Fig. 7a we show the to- 
tal emitted energy <§ [Eq. (35)] for three different values 
of tg; for the shortest pulses one can clearly observe an 
oscillating behavior of the total emitted energy as a func- 
tion of <3g/(27t); this behavior is less pronounced for the 
longest pulses. It is also worthwhile to notice that, for 
increasing values of Qg , the intensity can become so high 
that also for the shortest pulses the system is in any case 
completely destroyed by photoionization within the du- 
ration of the pulse itself. The increasing importance of 
photoionization implies a less remarkable difference in the 
time evolution of systems driven by pulses whose area is 
an odd or even multiple of n and, consequently, resonance 
fluorescence spectra characterized by a lower dependence 
upon the area of the pulse. 

In Fig. 76 we display tt-TaiS^x) I^i where S(u>x) is the 
central maximum value of the spectrum of resonance flu- 
orescence. The constant prefactor ttFa/2 allows us to 
compare the shape of the spectrum of resonance fluores- 
cence with that of a Lorentzian function of Auger de- 
cay width 7"a- If the only process involved was a decay 
causing a rate width P, then the spectrum of resonance 
fluorescence would be a Lorentzian function 



r/(2n) 



(38) 
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with peak value Lq = L(wx) and with total emitted en 
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Figure 8. Average resonance fluorescence spectrum for 500 
different realizations of the Gaussian driving pulses [Eq. (36)]. 
The mean duration of the pulses is tg = 7 fs, and the mean 
peak intensity is Ig = 7 x 10 17 W/cm 2 . Here tg and 7g 
are Gaussian random variables independently chosen for each 
realization with a variance equal to the 20% of the respective 
mean values. 



Figure 7. (Color online) (a) Total emitted energy S — 
J-^ S z (u), Q) du) and (b) peak value of the spectrum S(u)x) 
multiplied by 71-/^/2 as functions of the normalized pulse area 
Qg/(2-7t) [Eq. (37)] for tg = 2 fs (red, dashed line), tg = 5 fs 
(black, dotted line), and tg = 10 fs (green, solid line). 

ergy $l — f^°° L(ui) duj = ttFLq/2. By computing in 
Fig. 7b the quantity 7t/a<S'(wx)/2, we can relate it to 
the actual total emitted radiation of Fig. 7a and un- 
derstand the relative importance of Auger decay in re- 
lation to the other decay processes. By comparing the 
oscillating features in Fig. 7b with those of Fig. 7a. one 
notices that ttTaS^x) /2 approaches S only for short 
pulses satisfying Qq = 27r(n— 1/2). In these cases, as we 
have already discussed in Fig. 5a, the main term is repre- 
sented by post-x-ray-exposure Auger decay of the system. 
Nonetheless, because of the non-negligible role played 
by Rabi flopping, photoionization and spontaneous de- 
cay, one can notice in Fig. 7 a clear difference between 
irr A S{uj x )/2 and S. 

Figures 4 and 6 reveal that Rabi flopping produces a 
clear signature in the spectrum of resonance fluorescence 
of Gaussian pulses, which are becoming available by self- 
seeding at LCLS [58, 59]. However, since shot-to-shot 
variations in pulse intensity and duration are anticipated, 
we investigate how the spectrum of resonance fluores- 
cence is influenced by the presence of these fluctuations. 
For this purpose, we compute the energy spectrum of 
resonance fluorescence for a wide set of Gaussian pulses 
[Eq. (36)], by independently randomizing their duration 
and energy. The mean duration is chosen to be tq = 7 fs 
and the mean peak intensity is Iq = 7x 10 17 W/cm 2 , giv- 
ing a mean peak Rabi frequency of ss 6 eV. We compute 
the energy spectrum of resonance fluorescence for 500 



different realizations of the driving pulses. Thereby, the 
duration and the intensity are random variables whose 
probability distribution is Gaussian with a variance of 
20% of the mean value. The resulting average resonance 
fluorescence spectrum is shown in Fig. 8. It reveals that 
Rabi flopping is discernible even if the energy and dura- 
tion of the pulse vary appreciably from shot to shot. 



B. SASE x-ray pulses 

After investigating resonance fluorescence with laser- 
like regular Gaussian pulses, we turn to the presently 
available SASE pulses at LCLS. The SASE light is mod- 
eled with the partial coherence method (PCM) [93, 94], 
whose details are discussed in Appendix B. The SASE 
pulses have a central photon energy which is tuned to 
the transition energy of Ne + of 848 eV, with a band- 
width (FWHM of \£(uj)\ 2 ) of Aw SASE = 6 eV. The enve- 
lope function f(t) that we adopt [Eq. (B8)[ has FWHM 
duration r cnv = 6.5 fs. Further details are discussed in 
Appendix B. 

In Fig. 9 a we display the time-dependent Rabi fre- 
quency [Eq. (8)] J?r(£) = p£o(t) induced by the ampli- 
tude of a SASE pulse and in Fig. 9b the phase fx{t) 
of a SASE LCLS pulse obtained with the PCM method. 
The mean Rabi frequency and phase are also given. In 
Fig. 10a the time evolution of the population of the ex- 
cited state and the total population of the two-level sys- 
tem are plotted if the phase of the pulse is supposed 
constant, <px(t) = 0, and the spiky time-dependent Rabi 
frequency of Fig. 9a is used to integrate the EOMs. In 
Fig. 106 both the Rabi frequency and the phase of Fig. 9 
are used to integrate the EOMs. If the phase fluctua- 
tions are neglected, the decay of the system is slower; 
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Figure 9. (Color online) (a) The Rabi frequency J?r(/;) in- 
duced by the amplitude of a SASE pulse and (6) the phase 
¥>x(£) of the SASE pulse (red, solid lines) and their mean 
value (black, dashed lines). The mean pulse has a duration 
T cnv = 6.5 fs and a peak intensity / = 3.8 x 10 18 W/cm 2 . Its 
bandwidth is Acjsase = 6 eV. 



in both cases, due to the chaotic SASE pulse shape, the 
time evolution is very irregular. For the case displayed 
in Fig. 10a, though, one can see the presence of complete 
oscillations in p22(t), reaching its minimum at P22(t) = 
and its maximum when pn(t) = 0: this feature disap- 
pears when the phase fluctuations of Fig. 96 are taken 
into account. 

In contrast to the case of a Gaussian pulse, one can- 
not extract from the time evolution of the system any 
clear relation to the pulse area. Nevertheless, one ob- 
serves a relation between the Rabi frequency of Fig. 9a 
and the frequency with which the population of the ex- 
cited state P22(t) oscillates in Fig. 10. These oscillations, 
in fact, take place in a time interval which is shorter 
than the time characterizing the random fluctuations of 
Fig. 9. They are Rabi oscillations induced by the interac- 
tion with the intense driving field; as we show in Fig. 11, 
if a pulse of similar bandwidth but of far lower intensity 
is used to excite the system, the time evolution of the 
atomic system displays slower oscillations, whose mean 
frequency increases at increasing intensities. We further 
notice that, because of photoionization, the increase in 
the intensity of the driving field reduces the actual de- 
cay time of the system: this emerges by comparing the 
graphs displayed in Fig. 11 with that of Fig. 106. 

In Fig. 12 we display the resonance fluorescence spec- 
trum from SASE x rays. To observe Rabi flopping we 
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Figure 10. (Color online) Time evolution of a two-level sys- 
tem driven by a SASE pulse with the time-dependent Rabi 
frequency of Fig. 9a. The phase is assumed to be (a) constant, 
V?x(i) = 0, or (b) to be equal to the phase of Fig. 96. The 
red, dashed line shows the evolution of the total population 
Pn(i) + P22(t) [Eq. (9)]; the black, solid line represents the 
occupation of the excited state p22(t). 
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Figure 11. (Color online) Time evolution of a two-level system 
driven by SASE pulses generated with the PCM described in 
Appendix B. In both cases the pulses have a mean duration 
Tcnv = 6.5 fs and a bandwidth Aojsase = 6 eV. The peak 
intensity is (a) / = 3.8 x 10 15 W/cm 2 and (6) I = 8.8 X 
10 17 W/cm 2 . The red, dashed line shows the evolution of the 
total population pn(t) + p%a(t) [Eq. (9)]; the black, solid line 
represents the occupation of the excited state p22{t). 
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Figure 12. (Color online) Resonance fluorescence spectrum 
for SASE pulses. The black, dotted line shows the spectrum 
from a Gaussian pulse with FWHM duration tg = 6.5 fs and 
peak intensity Iq = 3.8 x 10 18 W/cm 2 . The red, solid line 
is the arithmetic mean over 1000 SASE pulses with average 
peak intensity Ig, FWHM duration tg, and a bandwidth of 
Acjsase = 6 eV. The green, dashed line is for the pulse in 
Fig. 9. 
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Figure 13. (Color online) Average resonance fluorescence 
spectrum over SASE pulses. The pulses have an average 
peak intensity of I = 1.6 x 10 18 W/cm 2 and bandwidth 
Acjsase = 6 eV. The red, solid line shows the average over 
SASE pulses with average FWHM duration of r en v = 6.5 fs; 
the black, dotted line is associated with pulses of average 
FWHM duration of r env = 2 fs. 



need the Rabi oscillations to occur within the coherence 
time of the pulse; for this reason, the intensity of the ex- 
ternal electric field is chosen such that the maximum Rabi 
frequency is larger than the bandwidth Awsase of the 
pulse itself. We look, in particular, at the emitted spec- 
trum by averaging over 1000 independent SASE pulses. 
The tails appearing in the spectrum of Fig. 12 are non- 
vanishing contributions at frequencies higher than the 
bandwidth of the pulse itself. These tails would not ap- 
pear if the field had equal bandwidth but lower intensity: 
they represent, therefore, a signature of the Rabi oscilla- 
tions described in Fig. 10. These tails are also in good 
agreement with the spectrum emitted when a Gaussian 
transform-limited pulse of identical intensity and time 
duration — but clearly with much lower bandwidth — is 
used to excite the system. If the phase of the SASE 
pulse remained constant and only its amplitude displayed 
chaotic fluctuations, then the spectrum emitted after one 
single pulse would be symmetric; furthermore, the aver- 
age spectrum would present a lower width, due to the 
absence of phase fluctuations. A clear observation of the 
tails of Fig. 12 and of the enlargement of the resonance 
fluorescence spectrum at increasing intensities might rep- 
resent a possible way to detect Rabi flopping also at 
present SASE facilities. 

Analogous conclusions had been drawn for the reso- 
nant Auger electron spectrum [51]: the width of the res- 
onant Auger electron line profile was expected to help in 
estimating the presence of Rabi oscillations in the system. 
Nonetheless, the very short coherence times at present 
XFEL facilities limited the actual experimental observ- 
ability of this effect at LCLS [52]. 

As a last point, we study the dependence of the reso- 
nance fluorescence spectrum on the duration of the SASE 



pulse. In Fig. 13 we plot the average spectrum emit- 
ted by Ne + cations when excited by an ultrashort pulse 
with peak intensity I = 1.6 x 10 18 W/cm 2 and band- 
width Awsase = 6 eV. The results are obtained by 
averaging over spectra resulting from SASE pulses, re- 
spectively, with a FWHM duration of r en v = 6.5 fs and 
of r env = 2 fs. It is worth noticing a remarkable differ- 
ence between different pulse durations. Naively, after the 
previous considerations, we would assume that the res- 
onance fluorescence peak has a FWHM associated with 
the large bandwidth of the pulse Awsase = 6 eV. For 
the shortest pulses, though, the resonance fluorescence 
spectrum exhibits a higher central peak whose width is 
clearly lower than Awsase- The explanation is based on 
the same arguments that we used to describe the spec- 
tra depicted in Fig. 6a, in which the post-x-ray-exposure 
decay results in a high Lorentzian peak of width given 
by the Auger decay width of the system. Analogously, 
for the ultrashort SASE pulses with r onv = 2 fs used in 
Fig. 13, the interaction with the pulse is shorter than the 
time needed by the system to completely decay; hence, at 
the end of the pulse, the probability of destruction of the 
system is about 90%. The Auger decay which follows 
the interaction with an ultrashort SASE pulse implies, 
therefore, the high central peak in the resonance fluores- 
cence spectrum shown in Fig. 13; for the same reason, 
its width is lower than the bandwidth of the pulse itself. 
A similar reduction of the FWHM was also observed in 
[52] in the Auger electron spectrum. In that case, by us- 
ing the same x-ray pulse to create Ne + ions and to drive 
the \s2p~ 1 — > Is -1 2p transition, the system could not 
completely Auger decay before the end of the pulse. The 
decay of the excited state with the natural decay time of 
the system turned out to dominate the observation. 
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IV. CONCLUSION 

In this paper we study theoretically the resonance flu- 
orescence for intense ultrashort x rays. The fundamental 
role played by resonance fluorescence for the study of the 
quantum properties of light renders it a cornerstone of 
x-ray quantum optics. It also gives an alternative and 
more easily available point of view on resonant Auger 
decay in ultraintense x rays, whose study, at present, is 
not yet fully conclusive [52]. Therefore, we investigate 
the nonlinear phenomena of Rabi flopping, i.e., repeated 
cycles of stimulated emission and absorption of photons 
induced by the interaction with the ultraintense pulses 
from XFELs, and its signature in the resonance fluores- 
cence spectrum at x-ray frequencies. 

We develop a two-level model of resonance fluores- 
cence whose time evolution is described by master equa- 
tions which include the coherent interaction of the system 
with the classical x-ray field. All processes that destroy 
the system, namely, Auger decay and photoionization, 
are fully taken into account. We use our model to de- 
scribe Ne + cations driven by an intense linearly polarized 
x-ray field tuned to the Is 2p _1 — » Is" 1 2p transition at 
848 eV; the transition is well isolated, i.e., separated by 
more than 70 natural linewidths from the lowest lying 
Rydberg excitation, Is — > 3p [52]. The intensity avail- 
able at present x-ray FELs such as LCLS is sufficiently 
high to induce Rabi flopping at frequencies that com- 
pete with the rate of destruction of the system. The 
two-level approximation allows us to investigate the res- 
onance fluorescence of photons associated with the tran- 
sition to the state with Ml = for two different sce- 
narios. First, we consider SASE radiation from present 
XFELs; second, we explore resonance fluorescence from 
coherent Gaussian pulses which are becoming available 
via the use of self-seeding techniques at fourth-generation 
x-ray sources. The measurement of the spectra predicted 
in this paper need to take advantage of the polarization 
properties of the emitted light. 

In the case of laserlike Gaussian pulses a clear signa- 
ture of Rabi flopping is predicted. We further show that 
the observation of Rabi flopping persists even when in- 
tensity and duration of the pulse vary appreciably from 
shot to shot. For SASE pulses, even though Rabi flop- 
ping does not manifest itself as clearly as in the previous 
case, we predict the appearance of tails in the spectrum 
that might represent a good signature of Rabi oscillations 
in the atomic system. These tails would not appear if the 
system was excited by a less intense pulse of equally large 
bandwidth. In the case of the resonant Auger spectrum, 
however, the presently large bandwidth at LCLS repre- 
sented a limit for the observation of analogous effects in 
the resonant Auger electron line [52] and the signature 
of Rabi flopping did not appear indistinguishably. Also 
in the case of resonance fluorescence the identification of 
Rabi flopping in the spectrum might be challenging. The 
amplitude of the aforementioned tails in which one is in- 
terested is predicted to be neither very high nor easily 



distinguishable. A much clearer signature is, however, 
identified for ions driven by Gaussian pulses, making the 
prospects with self-seeded LCLS very promising. 

The results which have been presented motivate fur- 
ther experimental investigation of resonance fluorescence 
at XFEL facilities. In particular, the method which has 
been discussed here is a good candidate for further stud- 
ies at hard x-ray frequencies. In the case of argon cations, 
for instance, the spectrum of resonance fluorescence, be- 
cause of the higher fluorescence yield compared with that 
of neon, is more intense than the one predicted in this pa- 
per. By including the radiative decay width in the EOMs 
(18), the model can be used to study the resonance flu- 
orescence spectrum of cations, e.g., argon, with higher 
fluorescence yield than the cations considered here. The 
basic features are discussed in this work and no qualita- 
tive differences are expected in heavier cations. Studies of 
resonance fluorescence in different atomic systems might 
require the use of a generalized formalism. The two-level 
approximation adopted in this paper, in fact, is not al- 
ways sufficient to properly describe the atomic transitions 
of interest: in these cases multilevel systems have to be 
considered. In addition, in order to take into account the 
atomic properties for different x-ray transitions, a con- 
siderable amount of new atomic data is necessary, such 
as, for example, Auger pathways and decay rates. A de- 
tailed theoretical study of the atomic properties of the 
system would have to be implemented, motivating fur- 
ther research in this field. 

Studies of resonance excitations followed by if-shell 
photoionization are receiving a lot of interest also for 
their potential applications in the biomedical sector [95- 
98]: Even though present facilities are not available 
yet for medical applications, studies of resonance fluo- 
rescence of If -shell transitions might also significantly 
contribute to the development of such applications of 
XFELs. 

In addition, resonance fluorescence plays a crucial role 
in the study of the nonclassical properties of light, such 
as photon antibunching [11, 12, 99, 100] and squeezing 
[101-104]. Our study opens the x-ray regime up for quan- 
tum optical effects which can be investigated by means 
of ultraintense pulses now available at XFELs. 

Finally, for nonstationary light, e.g., when the elec- 
tric field has a pulse-shaped envelope, the study of the 
time-dependent power spectrum [24, 77, 105-107], i.e., 
the time-dependent rate of detected photons, would allow 
one to investigate how the spectral properties of the fluo- 
rescent light evolve during the pulse. Even though such a 
power spectrum cannot be measured at present because 
of the ultrashort nature of XFEL pulses, with duration 
of the order of 10-100 fs, and because of the lack of suf- 
ficiently fast detectors, the study of the time-dependent 
power spectrum might provide better understanding and 
further knowledge of the interaction between matter and 
x rays. 
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Appendix A: Polarization effects and measurement 
geometry 



polarized along the axes e y and eg{9). Conversely, the 
photons spontaneously emitted to the state | lo) are ex- 
clusively polarized along the axis eg{9) and their inten- 
sity, varying in space as cos 2 6, is maximized for 6 = 0, n. 
For the same angles the intensity of the photons that are 
emitted in transitions to the two undriven states |1±) and 
which are linearly polarized along eg (9) vanishes. This 
motivates our choice throughout the paper of studying 
the spectrum of resonance fluorescence for 9 — 0, e r = e x 
and eg = e z . Polarization-dependent detection of the 
resonance fluorescence spectrum can take advantage of 
the properties just presented. 



We consider here the resonance fluorescence spectrum 
emitted by the two-level system displayed in Fig. 2 which 
is measured by rotating the detector around the y axis 
of Fig. 1, i.e. the spectrum at point r = re r (9), with 
e r (9) = cos9e x + sin9e z , where 9 is the angle between 
e r {9) and the x axis, lying in the x — z plane. We further 
introduce the vector eg(9) = — sm9e x + cos6>e z , which 
also lies in the x — z plane and is orthogonal to e r {ff)\ in 
this way, from (6) and (27), one has that E + (r e r (9), t) = 
E+{re r (9), t) e e {9) + E+ (r e r {9), t) e y , with 



E+(re r (9),t) 



suit 



cos0<7i o2 (*') 
<n +2 (i')- *i_a(f) 



V2 



(Al) 



and 



E+(re r (6),t) = -L ^1 (* 1+2 (f) + ^_ 2 (f)) , 

(A2) 

with t' = t — r/c. The autocorrelation func- 
tion (25) is G( l \h,t 2 ,re r {6)) = G { e l \tr,t 2 ,r e r (6)) + 
Gi 1) (t 1 ,t 2 ,re r (6)), with 

Gf ] {t u t 2 ,r e r {0)) = 2(r) [cos 2 9 (a 21o (t{) a lo2 (t' 2 )) 



((<T21_ &1_2(4)) + (*21 + (4) ^1+2(4))) 



and 



(A3) 



GV{t x ,t 2 ,re r {6)) = X(r) \ {{a 21 _ a x _ 2 (t' 2 )) 

+{^2i + (t' 1 )ai +2 (t 2 ))), 

(A4) 

where T(r) is defined in (33) and the application of the 
quantum regression theorem allows one to show that the 
crossterms (<X2i(*i) & 32 {t' 2 )), with i, j e {1+, 1-, lo}, i ^ 
j, vanish for any t\ and t' 2 . An analogous spatial depen- 
dence can be displayed also in the resonance fluorescence 
spectrum 5* (co, r e r {9)) — Sg(uj 7 r e r (9)) + S y (uj, re r {9)). 

In conclusion, we notice that the photons emitted in 
transitions to the two undriven states |l-j-) can both be 



Appendix B: Partial Coherence Method 

We use the partial coherence method (PCM) intro- 
duced in Ref. [93] to generate random realizations of the 
temporal shape of SASE XFEL pulses, whose knowledge 
is an important prerequisite for meaningful investigations 
of nonlinear x ray-matter interaction [51]. Those param- 
eters which can be measured at present XFELs, such as 
the average spectral intensity and the pulse duration, are 
taken into account as input parameters [94]. The PCM 
is used to generate non-transform-limited pulses, with a 
coherence time lower than the average FWHM duration 
of the pulse and with significant fluctuations in the pulse 
shape from shot to shot. 

The pulses are generated starting from their frequency 
representation £(w), whose amplitude is given by the av- 
erage spectral intensity of the pulse. If the phase of £(oj) 
was constant, then by Fourier transform one would ob- 
tain a transform-limited pulse. In order to generate a 
SASE pulse, we let the spectral phase vary in [— it, ir[. 

The PCM models the classical electric field £(t) 
[Eq. (1)]. We introduce the complex electric field [71] 
E ± (t) = \£ Q {t) eT'IwtWM and the complex field enve- 
lope^) = i£ We" iw(f) , such that £ (t) = £{t) c~ iuixt + 
£*(t)e iuxt . It follows that [71] 



\£(t)\ 2 = 2^(1)1" = 2\£(t)\' 



\£o(t)\ r < 



(Bl) 



We define the Fourier transform of £{t) as 

p + ca 

£{w)= £{t)c wt At = I^HIc-^ (B2) 



and from Parseval's theorem it follows that 



+oo 1 

2 j + _ 1 / I cf . .M2 



W d *=^/ |£M| 2 dw. (B3) 



Analogously one can define £(u) and E + (w) as the 
Fourier transforms of £(t) and E + {t) respectively. One 
notices that E + (uj) = £ (oj — uix) and therefore £(u>) = 
£ (u> — ujx)+£ (u) + wx), so that, from Parseval's theorem 
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(B3), one finds in agreement with (Bl) that 

1 f + OO 



1^)1 di = —/ 2\E(u>)\'6u>. (B4) 



The average spectral intensity of a SASE pulse is mod- 
eled here — close to measured spectral intensities — as a 
Gaussian function, so that 



|£H| 2 = |£ sp , max | 2 c-^ 2 / r22 ), 



(B5) 



whose FWHM is Aw S ase = 2/2Vln2. The FWHM dura- 
tion of the squared modulus of the inverse Fourier trans- 
form of \£ (oj)\ [71], which is here also a Gaussian function, 
is tsase = 41n2/AcjgASE- It follows that 



\£ HI 2 = |£„ x | 2 (c-^ x > 2 / fi2 J +e-t(^) 2 /^ 2 ]). 

(B6) 

The average spectral intensity, though, does not pro- 
vide any information about the spectral phase of the 
pulse. In analogy to the phase retrieval in x-ray crys- 
tallography [76], the knowledge of the spectral ampli- 
tude is not sufficient to completely determine the tempo- 
ral shape of the pulse via inverse Fourier transform. In 
the PCM approximate phase retrieval is achieved by as- 
suming initially a random frequency-dependent spectral 
phase varying in [— 7r, 7r[. 

We define a discrete spectral component of the elec- 
tric field £(u!i) = \£(u!i)\ e -1 ^ , with a sampling inter- 
val |wi-fi — coi\ <C AwgASE- 0i are random numbers in 
[— 7r, 7r[. The discrete inverse Fourier transform of £(u>i) 
provides the time-dependent discrete field lZ(tj). 

The complex function lZ{t), obtained by interpolating 
lZ(tj), spans an infinitely long interval in time because 
of the fluctuating <fii . To generate SASE pulses of finite 
duration, we multiply lZ(t) by a temporal filter function 
f(t). This function is non-zero only in a finite domain 
and the FWHM duration of \f{t)\ 2 is r onv . The finite du- 
ration of FEL pulses is determined by the electron bunch 
duration [108] and is usually measured. All together, we 
approximate the complex electric field by 



£(t) = -£ (t)e-^ = 1l(t)f(t). 



(B7) 



Along the way, we notice that the inverse Fourier trans- 
form of 7Z(t)f(t) is given by the convolution of the respec- 
tive inverse Fourier transforms £{lo) and f(uj). £ (to) has 



a random fluctuating phase <fi(uj), whereas |/(w)| 2 has 
a spectral FWHM Aw onv related to the inverse of r cnv . 
Hence, the spectral amplitude of a single pulse generated 
with the PCM also displays a spiky structure, where the 
average FWHM frequency of each spike is about Au; cn v 
[108]. In addition, since the average value of 4>(uj) is 0, 
the average spectral amplitude results from the convolu- 
tion of \£ (u)\ and /(w) and, because r onv 3> tsase 5 the 
width of f{ijj) is much narrower than the width of \£ (u))\. 
Consequently, the convolution 

I^-^OI/V) <l/«|£(«)|, 

-oo 

i.e., the multiplication by the envelope function fit) does 
not significantly affect the average spectral intensity of 
£(t). 

To generate SASE pulses for this paper (Fig. 9) we use 
the envelope function 



focas^t/T) if|,|<T/2 
\0 if|i >T/2 



with T = 7rr onv /(2 arccos \/l/2) and r cnv = 6.5 fs, de- 
fined as the FWHM duration of |/(t)| 2 [109]. The Fourier 
transform of f(t) is 



T sine 
/M = 77 fo 



2 J 



(B9) 



Then, Aw cnv « 2.41/r cnv is the FWHM of |/(w)| 2 . One 
notices that, for r onv = 6.5 fs and Awsase = 6 eV, one 
has Aw cnv = 0.24 cV < Aoj sase = 6 eV. 

Alternative approaches have also been developed and 
adopted, e.g., in [56, 110]. In these cases, the electric 
field is written as a Fourier series in time domain 



£{t) = ^2 a k cos(cjfet) + b k sin(wfct), 



(B10) 



k— — oo 



where the real coefficients and bk are independent 
zero-mean Gaussian random variables. Basically, this 
represents only a different description of £{lo) compared 
with the PCM. 
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